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Table 1: Run parameters and autocorrelation times for the topological charge for the simulated tree-level- 
improved Symanzik gauge ensembles. The lattice volume is kept hxed to L Ri 1.6 fm. The separation 
between the saved configurations of ^ 2 x T,„, {Q^) allows for the autocorrelation between the saved conhg- 
urations to be neglected. 


1. Introduction 

To determine potential channels for New Physics (NP) it is necessary to test and over-constrain 
the Standard Model (SM) by performing model independent theoretical predictions with high pre¬ 
cision. For this all systematic errors need to be controlled. Currently the heavy quark sector (charm 
and bottom) is rather unexplored compared to the light quark sector (see [|^) which makes it in¬ 
teresting for numerical simulations in Lattice Quantum Chromodynamics (LQCD). However, the 
simultaneous simulation of light and heavy (charm) quarks poses a difficulty. Namely, in order to 
control finite volume effects for the light quarks large lattice volumes are required, whilst to resolve 
the heavy quarks it is necessary to have fine lattice spacings a, which is very costly to achieve at 
the same time. The work presented in this talk investigates the feasibility of simulating charmed 
mesonic quantities within a Domain Wall formalism to lay the ground work for RBC/UKQCD’s 
ongoing 2+1/ charm program at the physical point whilst maintaining chiral symmetry and hence 
0{a) improvement |^. 

To this end we produced 4 tree-level improved Symanzik gauge ensembles with inverse lattice 
spacings between 2GeV and 5.6GeV in an approximately constant volume allowing for a 
controlled continuum limit. We stress that this quenched set-up is not intended to make physical 
predictions, but to have a feasible framework to test the applicability of Mdbius Domain Wall 
fermions to simulations of heavy quark physics. This is done by investigating the scaling behaviour 
in the continuum limit approach. To reduce the cost and since the simulated volume is comparably 
small, no light quarks are simulated but instead we consider q^-like D^-like and T]c-like quantities. 
Using the quenched approximation allows the exploration of very fine lattice spacings which are not 
otherwise affordable. Furthermore we developed the expertise and understanding of the parameter 
space which is needed for the dynamical simulations and are now set-up for the aforementioned 
ongoing 2+ 1/ measurements presented at this conference [^. A similar study comparing different 
valence discretisations for the heavy quarks was carried out on the same gauge ensembles and also 
presented at this conference |^]. 

2. Ensembles 

We perform all our measurements on four sets of ^(100) independent, tree-level-improved 
Symanzik gauge configurations 0, ||], where the inverse lattice spacing is varied in the range 2 — 
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Figure 1: Left; Determination of the lattice spacing with the Sommer scale parameter ro and flow parameter 
fo as a cross check of the determination of the lattice spacing. Right: log a as a function of j3 to determine 
finite size effects. The larger symbols correspond to the 4 ensembles used for measurements, the smaller to 
the auxiliary ensembles with volumes of 1.23 — 1.89fm. 


5.6 GeV and the physical volume is kept fixed to ~ 1.6fm (cf. table ||). For the generation of SU(3) 
gauge configurations, we use the heat-bath (HB) algorithm implemented in CHROMA [^] as well 
as additional code optimised for BG/Q. In Table |I| we give the parameters of the performed runs 
for the configuration generation and the autocorrelation times of the topological charge Q, 

Q= ( 2 . 1 ) 

where Ppv represents a clover definition of the field strength tensor. The topological charge was 
closely monitored to ensure ergodic sampling. The square topological charge is considered 
to be a slowest mode relevant for the dynamics of the simulated system and we use its auto¬ 
correlation time (found as described in [^]) to estimate the minimal separation between the saved 
configurations, so that they can be considered independent. 


2.1 Scale setting 


We use the Wilson flow to set the scale of our simulations, in particular, we use the wq - 
scale which was proposed in [10], From this we also take the physical scale = 0.1755(18)(04)fm, 
computed for 2 -|- 1/ simulations. The program package GLU^ used for the Wilson flow measure¬ 
ments was kindly provided by Jamie Hudspith. We additionally perform a cross check of the scale 
setting procedure by measuring the Sommer scale parameter vq and an additional flow parameter 
to 0 on the 3 coarsest ensembles. It is found that these agree well (cf. l.h.s of figure H). The 
one-loop beta function was used to estimate the parameters of the simulations. From the r.h.s. of 
figure H it can be seen that all lattices with volumes larger than ~ 1.2fm (i.e. all points other than 
the grey squares in figure [I]) are consistent with the fitted scaling curve. The universal scaling curve 


for all simulated volumes larger than 1.2fm confirms the observation of ref. [10| that the finite size 


'https://github.com/RJhucispith/GLU 
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Figure 2: Left: Behaviour of /ps y/nip^ for heavy-heavy pseudo scalar mesons as a function of the inverse 
pseudo scalar mass mps for different values of M 5 on the coarsest ensemble using Shamir DWFs with Ls = 
16. The data is normalised at mpg™ = l.SGeV (horizontal and vertical green dotted lines) to avoid finding 
a renormalisation constant. The vertical red line corresponds to the r\c mass. Right: Overlay of data from 
the two coarsest ensembles for 3 choices of M 5 . The dotted (dashed) curves depict data from the (second) 
coarsest ensemble. 

effects of the particular definition of the Wilson flow parameter are negligible in the volume range 
we are interested in. 

3. Parameter Choices 

In this study we aim to determine the range of parameter space in which Domain Wall fermions 
(DWF) can be used to reliably simulate heavy quarks with 0(a) improvement. The use of domain 
wall fermions introduces two new parameters which can be chosen freely: The Domain Wall height 
M 5 and the extent of the fifth dimension L^. In this study we use two different versions of the 
Domain Wall formalism both of which are used in the RBC/UKQCD 2 + 1/ simulations, namely 
the ’standard’ Shamir Domain Wall fermions [ |I^ |I^ and the more recently developed Mobius 
Domain Wall fermions [ 0 - 111 . Mobius Domain Wall fermions are a rescaled version of the 
Shamir case, allowing to simulate the same physics with half the extent in the 5th dimension and 
therefore decreasing the computational cost. In particular choosing _ ^Shamir ^2 is expected 

to reproduce the same results. 

When naively increasing the bare quark mass, unphysical behaviour is found for > 0.4. 

This can be seen from the unexpected bending down of /psy^mps as heavier pseudo scalar masses 
are approached (compare left hand panel of figure The right hand panel of figure ^ suggests that 
this behaviour is very sensitive to the choice of M 5 . Additional to the pseudo scalar mass for which 
this bending sets in, the size of the cut-off effects also depends strongly on the choice of M 5 . We 
found that the optimal choice of the Domain Wall height is given by M 5 = 1.6, minimising cut-off 
effects at the same time as postponing the unphysical behaviour to very close to charm even on the 
coarsest ensemble. In the limit Ls Domain Wall fermions maintain exact chiral symmetry on 
the lattice. However for finite Ls there is a residual chiral symmetry breaking that can be quantified 
by the residual mass /Wres which can be defined from the axial Ward Identity for the Domain Wall 
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Figure 3; Left; Behaviour of the effective residual mass as a function of time, m®® is found by evaluating 
(3.2) on every time slice. Right; The residual mass taken at the centre of the lattice as a function of the bare 
quark mass. Both panels are for the coarsest ensemble with Shamir Domain Wall fermions and Ls = 16. 


fermions 

{d^^^{x)PiO)) = 2m {P{x)P{0)) +2 {j 5 ,ix)PiO)) , (3.1) 


so that 


_ Lv/5,(x)P(0) 
LsPi^o) • 


(3.2) 


Here is the conserved axial current, P is the pseudo scalar density and J^g is a pseudo scalar 
current that mediates between the boundaries of the 5th dimensions and the centre of the 5th di¬ 
mension. 

To investigate the above mentioned unphysical behaviour the residual mass was closely moni¬ 
tored (compare figure H) and from this it can be seen that for amg > 0.4 the residual mass does not 
remain well controlled. For this reason we restricted the following Mobius Domain Wall fermion 
simulations to amg < 0.4 and effectively increased Ls by choosing it as = 12 (corresponding 

to = 24) as opposed to = 16 as in figures ^ and ||. 


4. Analysis and Results 

The physical quantities of which we investigate the continuum scaling are the decay constants 
of strange-strange (rjs), heavy-strange (D^-like) and heavy-heavy (T]c-like) pseudo scalar mesons. 
Our strategy is as follows; 

On each ensemble we simulate two closely spaced strange quark masses and various ’heavy’ 
quark masses with amg < 0.4 and compute pseudo scalar masses and matrix elements. Ds and 
rjs masses are interpolated to the physical strange quark mass, by matching the T]^ mass to the 
published value [|T7||. To be able to take a continuum limit at the same physical point the decay 
constants for rjc and are now interpolated to common reference masses mf^ on each ensem¬ 
ble. Since we are not interested in physical predictions in the quenched approximation we do not 
calculate renormalisation constants for these quantities, but instead take ratios at a reference mass 
to cancel the renormalisation constant. This reference mass is chosen as ntp”™ = 1 .OGeV. In sight 
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Figure 4: Continuum limit of the ratio of decay constants at different pseudo scalar masses for Ds (left) and 
Tjc (right) using Mobius Domain Wall Fermions. The data is normalised at mj)™™ = l.OGeV. We compare 
two different polynomial fit ansatze, namely linear (dashed lines, square-symbols) and quadratic (dotted 
lines, round symbols) polynomials in a^. For the heaviest two reference masses the coarsest ensemble does 
not allow a simulation with airiq < 0.4, so the continuum limit is taken from the three finer ensembles only. 

of further investigations towards the feasibility of B-physics and contact to HQET we will consider 
the expression fpsy/nips which approaches a constant in the static limit. The resulting continuum 
limit is presented in figure From the two panels in figure ^ it is clear that the continuum approach 
with the chosen parameters is very flat and well described by 0{a^) scaling. This is further encour¬ 
aged by the fact that the residual mass remained well behaved and under control for the simulated 
data points, confirming fhe 0(a) improvemenf. 

5. Conclusion and Outlook 

From fhese resulfs we conclude fhaf Mobius Domain Wall fermions are a suifable discrefisafion 
for charm physics. We expecf fhe qualifafive behaviour of our findings fo remain unchanged when 
going beyond fhe quenched approximafion allowing us fo dynamically simulafe QCD wifh 2 + 
If whilsf keeping discrefisafion errors under confrol. This is very encouraging for our ongoing 
dynamical efforls (firsl resulfs were presenfed in Jiiffner’s falk af fhis conference [§]). 

Finally fhis quenched dafa will be used fo fesf fhe feasibilily of doing B-physics wifh Domain 
Wall fermions, by combining if wifh resulfs in fhe sfafic limif as well as using fhe rafio mefhod [p^. 
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